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ABSTRACT 

Non-perturbative studies of quantum gravity have recently suggested the possibility that the 
strength of gravitational interactions might slowly increase with distance. Here a set of generally 

•rH , 

^ ■ covariant effective field equations are proposed, which are intended to incorporate the gravitational, 

H ; 

^ vacuum-polarization induced, running of Newton's constant G. One attractive feature of this 

approach is that, from an underlying quantum gravity perspective, the resulting long distance 
(or large time) effective gravitational action inherits only one adjustable parameter ^, having the 
units of a length, arising from dimensional transmutation in the gravitational sector. Assuming 
the above scenario to be correct, some simple predictions for the long distance corrections to 
the classical standard model Robertson- Walker metric are worked out in detail, with the results 
formulated as much as possible in a model-independent framework. It is found that the theory, 
even in the limit of vanishing renormalized cosmological constant, generally predicts an accelerated 
power-law expansion at later times t ~ ^ ~ 1/H. 
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1 Introduction 



Non-perturbative studies of quantum gravity have recently suggested the possibihty that gravi- 
tational couplings might be weakly scale dependent due to nontrivial renormalization group effects. 
This would introduce a new gravitational scale, unrelated to Newton's constant, required in or- 
der to parametrize the gravitational running in the infrared region. If one is willing to accept 
such a scenario, then it seems difficult to find a compelling theoretical argument for why the 
non-perturbative scale entering the coupling evolution equations should be very small, comparable 
to the Planck length. One possibility put forward recently is that the relevant non-perturbative 
scale is related to the curvature and therefore macroscopic in size, which could have observable 
consequences. One key ingredient in this argument is the relationship, to some extent supported 
by Euclidean lattice results combined with renormalization group arguments, between the scaling 
violation parameter and the scale of the average curvature. Irrespective of the specific details of a 
gravitational theory at very short distances, such results would bring gravitation more in line with 
the rest of the Standard Model, where all gauge couplings are in fact known to run. 

In this paper we investigate the effects of a running gravitational coupling G at large distances, 
with as few assumptions as possible about the ultimate behavior of the theory at extremely short 
distances, where several possible scenarios include a string cutoff at length scales A5 = (27rQ;')^/^ [1], 
the appearance of higher derivative terms (either as direct contributions or as radiative corrections), 
or perhaps a - somewhat less appealing - explicit ultraviolet cutoff at the Planck scale. The 
running of the gravitational coupling will generally be assumed to be driven by graviton vacuum 
polarization effects, which produce an anti-screening effect some distance away from the primary 
source, and therefore tend to increase the strength of the gravitational coupling. The above scenario 
is quite different from what one would expect for example in supergravity theories, where significant 
cancellations arise in perturbation theory between graviton and matter loops [2], and in contrast 
to ordinary gravity where in weak field perturbation theory L loops contribute L + 1 powers of the 
curvature tensor to the effective action [3]. Instead, the running of Newton's constant is thought 
to arise due to the presence of a non-trivial, genuinely non-perturbative, ultraviolet fixed point 
[4, 5, 6] (a phase transition in statistical mechanics parlance [4]). 

In this paper a power law (as opposed to a logarithmic) running of G will be implemented 
via manifestly covariant nonlocal terms in the effective gravitational action and field equations. It 
ultimately will involve the inverse of the covariant d'Alembertian raised to some fractional power 
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l/2z^, which in the framework of the present paper remains largely unspecified, although non- 
perturbative models for quantum gravity have recently put forward some rather specific predictions. 

Let us recall here, to provide some degree of motivation, the recent discussions of [7, 8] as 
a possible theoretical framework for the running of Newton's G. The above results suggest that 
the gravitational constant G cannot be regarded a constant as in the classical theory, but instead 
changes slowly with scale due to the presence of weak gravitational vacuum polarization effects, in 
a way described by 

G(r) = G(0) [ 1 + [r/if'' + 0{{T/ifh ] (l-l) 

The exponent generally related to the derivative of the beta function for pure gravity evaluated 
at the non-trivial ultraviolet fixed point via the relation (3'{Gc) = —l/v, is supposed to universally 
characterize the long-distance properties of quantum gravitation, and is therefore expected to be 
independent of the specifics related to the nature of the ultraviolet regulator, or other detailed 
short distance features of the theory. ^ 

Recent estimates for the value of the universal scaling dimension = —(3'{Gc) derived from 
non-perturbative studies of gravity vary from u"^ 3.0 [8] in the Euclidean Regge lattice case, to 
v"^ 3.8 in the 2 + e expansion [10] about two dimensions carried to two loops [11, 12, 13, 14], 
to ly-^ « 2.7 [15] and i^-^ « 1.7 [16] in an approximate renormalization group treatment a 
la Wilson based on an Einstein-Hilbert truncation, with some significant uncertainties in all three 
approaches. More details, as well as a systematic comparison of the various methods and estimates, 
can be found in [17], where we argued, based on geometric arguments, in favor of the exact value of 
the exponent = 1/3 for pure gravity in four dimensions, and 0{l/{d— 1)) for large d. It is perhaps 
a testament to how far these calculations have progressed that actual numbers have emerged which 
can meaningfully be compared between different (lattice and continuum) approaches. It should 
also be noted that, from a quantum gravity perspective, there are really no adjustable parameters 
in Eq. (1.1), except for the new non-perturbative curvature scale ^: both and v are in principle 
finite and calculable numbers. 

The mass scale m = in Eq. (1.1) is supposed to determine the magnitude of quantum 
deviations from the classical theory, and separates the short distance, ultraviolet regime with char- 
acteristic momentum scale /j, <^ m, where non-perturbative quantum corrections are negligible, 
from the long distance regime where quantum corrections become significant. The magnitude of ^ 
itself involves, in a rather non-trivial way, the dimensionless bare coupling G, the fixed point value 

^Already in ordinary Einstein gravity one finds for very sfiort distances r Ip corrections to the static potential, 
wliicli can be computed perturbatively [9]. In general for such short distances string corrections and/or higher 
derivative terms should be considered as well. 
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Gc and the ultraviolet cutoff A, 



— 1 



oc A exp 
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(1.2) 



Ultimately to make progress and determine tfie actual pfiysical value for the non-perturbative scale 
^ some physical input is needed, as the underlying theory cannot fix it (the ratio of the physical 
Newton's constant to can be as small as one desires, provided the bare coupling G is very close 
to its fixed point value Gc)- It seems natural to identify 1/^^ with either some very large average 
spatial curvature scale, or perhaps more appropriately with the Hubble constant (as measured 
today) determining the macroscopic expansion rate of the universe via the correspondence 



Let us briefly digress here, and recall that in non-Abelian SU{N) gauge theories a similar set 
of results is known to hold for the renormalization-group induced running of the gauge coupling g, 
so it will be instructive to draw further on the analogy with QCD, and non-Abelian gauge theories 
in general. Of course one crucial difference between gravity and ordinary gauge theories lies in the 
fact that, in the latter case, the evolution of the coupling constant can be systematically computed 
in perturbation theory due to asymptotic freedom, a statement which is known to reflect the fact 
that such theories become non-interacting at short distances, up to logarithmic corrections, making 
perturbation theory consistently applicable. It is well known that for weak enough gauge coupling 
in SU (N) gauge theories one has 



with /?o the coefficient of the lowest order term in the beta function, /i = 1/r an arbitrary momentum 
scale, Ajg^ ~ 220MeV a non-perturbative scale parameter, and the dots denoting higher loop 
effects. Instead of the Aj^g parameter one could just as well use some other physical scale, such as 
the inverse of the gauge correlation length, mQ++ = where the 0^^ denotes the lowest glueball 
state (the Slavnov-Taylor identities prevent of course the gluon from acquiring a mass to any order 
in perturbation theory). For the purpose of comparing to gravity, one should perhaps emphasize 
that confining non-abelian gauge theories such as QCD do not, and cannot, directly determine the 
scale which needs to be ultimately fixed by experiment from say a direct measurement of 

''a possible scenario is one in which = Hoc = limt^oo H{t) = \/0.a Ha with = = where A is the 

observed cosmological constant, and for which the horizon radius is -Roo = . 



(1.3) 



in a system of units for which the speed of light equals one. ^ 




(1.4) 
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the size of scaling violations. Its magnitude involves in a non-trivial way the bare gauge coupling 
g and the ultraviolet cutoff A, 

cx A exp J' (1.5) 

which is very much analogous to Eq. (1.1). The correspondence with QCD and nonabelian gauge 
theories would therefore suggest ^ Aj^g, with the gravitational ^ a new non-perturbative scale, 
ultimately also to be determined from experiment. 

Although not always necessarily advantageous (most perturbative calculations, being based on 
Feynman diagrams, are eventually done in momentum space and do not seem to benefit significantly 
from this approach), the running of the gauge coupling g can be re- formulated in terms of an effective 
action, involving the d'Alembertian acting on functions of the field strength. One sets 

with 2/3o = (lliV — 2n/)/(247r2) for non-abelian SU{N) gauge theories with Uf massless fermion 
flavors, and with the log of the d'Alembertian □ suitably defined, for example, via 

logf^) = rdm^l^ — ^ (1-7) 

V/"^/ Jo I 771,2 ^2 _,_ □ j V / 

leading to a one-loop corrected effective action of the form 

leff = i / dxF^,{x) (^^ + (3o log(;^) + •••) F'^'i^) (1-8) 

with jj, an appropriately chosen mass scale [18]. 

In the gravitational case the corrections described by Eq. (1.1) have a more complicated struc- 
ture, and in particular are no longer logarithmic. But they can be viewed for example as arising 
from a resummation of an infinite number of loop logarithms, as in the expansion 

j:i{-iT(^o,eor-{^)' (1.9, 

n=0 ^ 

In the next section we shall describe how the renormalization group induced running of the gravi- 
tational constant can be implemented in a simple way via a non-local set of manifestly covariant 
correction terms arising in the effective, long distance gravitational field equations. These effective 
equations can then be used as a basis for a systematic discussion of various quantum corrections 
to the standard solutions of the classical field equations. 
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2 Effective Gravitational Action and Effective Field Equations 



In general terms, a quantum-mechanical running of the gravitational coupling implies the re- 
placement 

G - G{r) (2.1) 

in classical physical observables. This is easier said than done, as in gravity the r in the running 
coupling G(r) is coordinate dependent, and as such can lead to considerable ambiguities regarding 
the interpretation of exactly which distance r is involved. A more satisfactory approach would 
replace G{r) in the gravitational action 

with a manifestly covariant object, intended to correctly represent an invariant distance, and in- 
corporating the running of G as expressed in Eq. (1.1), 

with the covariant d'Alembertian operator □ defined through an appropriate combination of co- 
variant derivatives 

□ = (2.4) 

Multiplication by the coordinate r gets therefore replaced by the action of □ , whose Green's function 
in D space-time dimensions is known to behave as 

1 1 

< x\-\y > 6{r - d{x,y\g)) ~ (2.5) 

Here d would be the distance along a minimal path z^{t) connecting the points x and y in a fixed 
background geometry characterized by the metric g^i,, and given by 

dix,y\g) = r^'\r^g^^[z)^^ . (2.6) 

Jt(x) 

As a result l/D can be envisioned as a coordinate independent way of defining consistently what 
is meant by r in the running of G(r), 

G(r) ^ G(n) (2.7) 

The above prescription has in fact been used successfully for some time to systematically incorporate 
the effects of radiative corrections in an effective action formalism [19, 20, 21]. It should be noted 
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that the coefficient in Eq. (1.1) is expected to be a calculable number of order one, but not 
necessarily the same as the coefficient cn, as r and l/y^ are clearly rather different entities to 
begin with. ^ 

One should recall here that in general the form of the covariant d'Alembertian operator □ 
depends on the specific tensor nature of the object it is acting on, 

° T'^'^-^S... = 9^"^,^ (v. T'^'^-^s..) (2.8) 
Thus on scalar functions one obtains the fairly simple result 

US{x) = J-d^g^'^'^d^Six) (2.9) 

whereas on second rank tensors one has the significantly more complicated expression DTq^ = 
g^'^'^ vTaji) ■ Furthermore one should recognize that the form for the effective gravitational 
action of Eq. (2.3) is possibly not unique. A more integration-by-parts symmetric expression would 
be, for example, 

' = 16^ /'"^ ^ - {-^T^ + • • •) ^ ^'-^^^ 

In general the covariant operator appearing in the above expression, namely 

^(°) = CO (^j (2.11) 

has to be suitably defined by analytic continuation from positive integer powers. The latter can be 
done either by computing for positive integer n and then analytically continuing to n — ^ —1/21^, 
or alternatively by making use of the identity 

1 (-1)' 



/ dss''-^exp(isn) (2.12) 
io 



r(n) 

and subsequent use of the Schwinger-DeWitt representation for the kernel exp(z s □) of the massless 
operator □. Within the limited scope of this paper, we will be satisfied with computing the effects 

of positive integer powers n of the covariant d'Alembertian □, and then analytically continue the 
answer to fractional n = —l/2v. In the following the above analytic continuation from positive 

integer n will always be understood. ^ 

^In the lattice theory cj was originally estimated from the invariant curvature correlations at around Cj « 0.01, 
while more recently it was estimated at cj w 0.06 from the correlation of Wilson lines [17]. It is important to note 
that while the exponent v is universal, in general depends on the specific choice of regularization scheme (i.e. 
lattice regularization versus dimensional regularization or momentum subtraction scheme). 

®We notice in passing that in this approach it is not obvious how to formulate a running cosmological constant, 
as the d'Alembertian □ in A(r) / dx^Jg — + A J dx^Jg (l — cu (l/C^^)^) has no function of the metric left to act 
on [22]. This situation is not entirely surprising as, lacking derivatives, the effect of the A term is just to control the 
overall scale. In pure lattice gravity the bare A is trivially scaled out and does not run [8, 17]. In this scenario the 
physical long distance cosmological constant ~ 1/C^, being related to an average curvature, is considered a physical 
quantity to be kept fixed as the gravitational coupling G(r) slowly evolves with scale. 
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It should be stressed here that the action in Eq. (2.10) should be treated as a classical effective 
action, with dominant radiative corrections at short distances (r ^ ^) already automatically built 
in, and for which a restriction to generally smooth field configrirations does make some sense. In 
particular one would expect that in most instances it should be possible, as well as meaningful, to 
neglect terms involving large numbers of derivatives of the metric in order to compute the effects 
of the new contributions appearing in the effective action. ^ 

A number of useful results can already be obtained from the form of the effective action in 
Eq. (2.10). In particular, once a specific metric is chosen, the running of G can be readily expressed 
in terms of the coordinates appropriate for that metric. Later in this work we will illustrate 
extensively this statement for the specific, and physically relevant, case of the Robertson- Walker 
(RW) metric. 

The next major step involves a derivation of the effective field equations, incorporating the 
running of G. As will be shown below this is not entirely straightforward, as the variation of 
the non-local effective action is complicated by the presence of a differential operator raised to a 
fractional power, acting on what are rather complicated functions of the metric to begin with. We 
shall therefore postpone a discussion of this aspect to the Appendix, which focuses on this specific 
topic. 

Had one not considered the action of Eq. (2.10) as a starting point for constructing the effective 
theory, one would naturally be led (following Eq. (2.7)) to consider the following effective field 
equations 

R^.u - \g^..R + ^9^.u = 8nG (1 + A{n)) T^, (2.13) 

the argument again being the replacement G{r) — G {1 + A{n)) involving the invariant object 
□ . Here, following common notation, A is the scaled cosmological constant, not to be confused 
with the ultraviolet cutoff. Being manifestly covariant, these expressions at least satisfy some of 
the requirements for a set of consistent field equations incorporating the running of G. The above 
effective field equation can then be easily re-cast in a form similar to the classical field equations 

Ri,u -\g^^vR + ^g^v = SttGT^, (2.14) 

with T/^i^ = (1 -|- A{\J)) T^i, defined as an effective, or gravitationally dressed, energy-momentum 
tensor. Just like the ordinary Einstein gravity case, in general T^i^ might not be covariantly con- 
served a priori, T^iy 7^ 0, but ultimately the consistency of the effective field equations demands 

'^Dominant contributions to the original Feynman path integral for the underlying quantum gravity theory are, 

on the other hand, presumably nowhere difFerentiable, the smooth configurations having ultimately zero measure in 
the gravitational functional integral [23]. Furthermore, issues related to causality, unitarity and positivity are better 
referred to the original, local microscopic action, which presumably shares all of these properties. 
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that it be exactly conserved in consideration of the Bianchi identity satisfied by the Riemann tensor. 
The ensuing new covariant conservation law 

V^f^. = [{1 + A{n)) T^,] = (2.15) 

can be then be viewed as a constraint on Tfj^i, (or T^u) which, for example, in the specific case of 
a perfect fluid, will imply again a definite relationship between the density p{t), the pressure p{t) 
and the RW scale factor R{t), just as it does in the standard case. 

This point is sufficiently important that we wish to elaborate on it further. In ordinary Einstein 
gravity the energy momentum tensor is defined via the variation of the matter action 

6lM = lj dx^Sg^^T^'' (2.16) 

But when the above arbitrary variation Sg/^i, is taken to be a gauge variation, 

= di,€^ + g\u 9^e^ + gfj,u (2.17) 

integration by parts in Eq. (2.16) immediately yields the covariant conservation law V^Tjj^i, = 0, 
as a direct consequence of the gauge invariance of the matter action. 

On the other hand, in the modified field equations of Eq. (2.13), the object which will be 
required to be conserved by the consistency of the field equations is the gravitationally dressed 
energy momentum tensor, namely (1 + ^4(0)) Tf^u, and not the original bare T^i, itself. Referring 
therefore to the original as "the energy momentum tensor" would appear to be improper, since, 
for the consistency of the effective field equations of Eq. (2.13), the latter is no longer required to 
be covariantly conserved. ^ In a sense, the effective field equations of Eq. (2.13) can be seen simply 
as a consequence of having changed the expression in Eq. (2.16) to 

5Im = \j dx^5g^, (l + ^(a)) r'^' (2.18) 

Let us make a few more comments regarding the above effective field equations, in which we will set 
the cosmological constant A = from now on. One simple observation is that the trace equation 
only involves the (simpler) scalar d'Alembertian, acting on the trace of the energy-momentum 

tensor 

R = SttG (l + ^(n)) (2.19) 

*This can be illustrated further by the specific case of the perfect fluid, for which the energy momentum tensor is 
usually written as = {p{t) + p{t)) u^u^ + g^^u pit). In general its covariant divergence is not zero, but consistency 
of the Einstein field equations demands V^T^v = 0, which for the RW metric forces a definite relationship between 
R{t), pit) and pit), namely pit) + 3 (p(t) [Rit)/Rit)^ = 0, irrespective of the equation of state relating p to 

p. In the effective field equations of Eq. (2.13) the perfect fluid form for T^j^ can still be used (at it still satisfies all 
the original symmetry requirements), but the covariant conservation law has the new form displayed in Eq. (2.15), 
which imposes a new constraint on the scale factor Ji(i), as well as on the underlying p(i) and pit). 
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Furthermore, to the order one is working here, the above effective field equations should be equiv- 
alent to 

(i - A{a) + 0{A{Df)) (r^, -ka^uR) = SttGT^, (2.20) 

where the running of G has been moved over to the "gravitational" side. Indeed it has recently 
been claimed [22] that equations similar to the above effective field equations (at least for positive 
integer power n, including the classical case tt, = 0) can be derived from a nonlocal extension of the 
Einstein-Hilbert action. In the classical case {A{0) = 0) one writes a new non-local action 

^ = 16^ / '^'^^ (^'^ ~ ^ ^) 5^'^'^ ^^-^^^ 
whose variation, it is argued, gives the correct field equations up to curvature squared terms 

V9 {R^'' -\g^^R + (rI,)) = (2.22) 
For non-vanishing A{n) the above construction can then be generalized to 

^ = 16^ / (^'^ ~ ^ ^) " + 0{A{nf)) ^R,, (2.23) 
whose variation can now be shown to give 

V5 (l - AiD) + OiAiDf)) (r^, -Ig^^R^ + O (rI,) = (2.24) 

and which would coincide with the previously proposed effective field equations, again up to higher 
order curvature terms. 

3 Covariant d'Alembertian on Scalar Functions 



As a first step in solving the new set of effective field equations, consider first the trace of the 
field equation in Eq. (2.19), written as 

(l - A{U) + 0{A{uf)^ R = %itGT/ (3.1) 

where R is the scalar curvature. Here we have made the choice to move the operator A{U) over 
on the gravitational side, so that it now acts on functions of the metric only, using the binomial 
expansion of 1/(1 + A{U)). A discussion of the full tensor equations and their added complexity 
will be postponed to the next section. To proceed further, one needs to compute the effect of A{U) 
on the scalar curvature. The d'Alembertian operator acting on scalar functions S{x) is given by 

^d^g^^-^d,S{x) (3.2) 
V 9 
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and for the RW metric, acting on functions of t only, one obtains a fairly simple result in terms of 
the scale factor R{t) 

As a next step one computes the action of □ on the scalar curvature R, which gives 

-6 [-2kR{t) -5R^{t) R{t) + R{t) R?{t) + 3R{t) R{t) R^^\t) + R^{t) /R^{t) (3.4) 

and then on i? which gives 

6 (-6 k R?{t) R{t) - 15 R!^{t) R{t) + 6k R{t) R'^{t) + 45 R{t) R^{t) R'^{t) - 
12 i?2 (t) i?3 {t)+6k R{t) R{t) (t) + 15 R{t) i?3 (t) (t) - 41 i?^ (t) R{t) R{t) (t) + 

5 R^ (t) ^{t)-2kR^{t) W {t)-9R^ (t) R^ (t) (t) + 7R^t) R{t) ^ (t) + 

4 R^t) R{t) (t) + i?(6) {t))/R{tf (3.5) 

etc. It should already become clear at this point that the computed expressions are rapidly becom- 
ing quite complicated. Nevertheless some of the higher order terms can, for example, be interpreted 
as higher derivative curvature contributions, since for Riemann squared, Ricci squared and scalar 
curvature squared, one has respectively 

Rf^uXa R'"'^" = 12 (k^ + 2kR\t) + R\t) + RHt)R\t)) /R{tf (3.6) 

R^^Rf''' = 12(A;2 + 2kR'^{t) + kR{t)R{t) + i?^(t) + R^{t)R^{t) + R{t)R^{t)R{t)) /Rit)"^ 

(3.7) 

R^ = 36 (A; + R^{t) + R{t)R{t)f /R{tf (3.8) 

with 

R-iivXa R^^^" — g Rfj^y R^^ — ^ i?^ = (3-9) 

for arbitrary scale factor R{t). But in the following we will just simply set R{t) = Rot°', in which 
case 



l^^XaR""^'' = ^ — (3.10) 



_ 12a^(2«^-2a + l) 

= (3.12) 

and for the scalar curvature (here allowing for k ^ 0, see Eq. (A.9) in Appendix A) 

k ^ a (—1 + 2 Q 
•-0 
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Acting with on the above scalar curvature now gives for k = 

6a {-l + 2a)t~'^ 



(3.14) 



36 (-1 + a) a (-l + 2a)r^ (3.15) 

144 (-1 + a) a (-1 + 2 a) (-5 + 3 q) (3.16) 

864 (-1 + a) a {-1 + 2 a) (-7 + 3 a) {-5 + 3a)r^ (3.17) 
for n = 0, 1, 2, and 3 respectively, and therefore for arbitrary power n 

c„6a (-l + 2a)r2-2" (3.18) 

with the coefficient given by 



r(^-n) ^ ^ 



^^ r(n + i)r(g^) 

Here use has been made of the relationship 

to analytically continue the above expressions to negative fractional n [24] . For n = —l/2v the 
correction on the scalar curvature term R is therefore of the form 

[l-c^it/O^/") ■ 6a {-l + 2a)t-^ (3.21) 

with 

^_i r(i - ^)r(^%^) 

= 2 - T^,3a-i , TV" (3-22) 
In particular for a = 2/3 (the classical value for a pressureless perfect fluid) and = 1/3 one has 

_ p-3 r(-^)r(^) _ TT 

- ^ p^2) - 4 ^^■''^^ 

whereas, for example, for a = 1/2 and u = 1/3 one obtains c^, = — \/7rr(|)/r(|). Putting 
everything together, one then obtains for the trace part of the effective field equations 

(l - CO c. (^'''^ + O (m?/^)^ ^-i2a-l) ^ 3^^^^^^ (3 2^^ 

The new term can now be moved back over to the matter side (since the correction is assumed to 
be small), in accordance with the structure of the original effective field equations Eqs .(2.13) and 
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(2.19), and thus avoids the problem of having to deal with the binomial expansion of 1/(1 + A{0)). 
One then has 

'-^S^p^ = 8.G (l + coc. (0''^ + O((t/0^/^)) Pit) (3.25) 

which is the RW metric form of Eq. (2.19). If one assumes for the matter density p{t) ^ p^t^ , then 
matching powers when the new term starts to take over at larger distances gives the first result 

/3 = -2 - \/v (3.26) 

Thus the density decreases faster in time than the classical value (/3 = —2) would indicate. The 
expansion appears therefore to be accelerating, but before reaching such a conclusion one needs to 
determine the time dependence of the scale factor R(t) (or a) as well. 

One might be troubled by the fact that some of the Gamma functions appearing in the expression 
for Cu can diverge for specific choices of i^, e.g. when v = l/2(n + 1) as in Eq. (3.22) for n integer. 
But further thought reveals that this is not necessarily a concern here, as the coefficient Cj, actually 
has to be divided out and then multiplied by (which, as discussed in the introduction and in 
[17], is expected to be a number of order one) to get the correct magnitude for the correction. One 
has therefore 

cp = (3.27) 

so that the correction eventually ends up as (1 + c^it/^)^/^), as it should, in accordance with 
Eq. (1.1) for G(r) (the "t" here is like "r" there). 

Having completed the calculation of the quantum correction term acting on the scalar curvature, 
as in Eq. (3.1), one can alternatively pursue the following exercise in order to check the overall 
consistency of the approach. Consider acting on T^^ = —pit) instead, as in the trace of the 
effective field equation Eq. (2.19) 

R = 8ttG (l + ^(n)) T^f" (3.28) 

for A = and pit) = 0. For pit) = po t^ and Kit) = Rof^ one finds in this case 

r(i. _|_ i^rf /3+3"+l \ 

□"(-/9(t)) ^ 4"(-l)"+^ . ' \J^/ pot^-^"" (3.29) 

^ ^ ^ r(f + 1 - n)r(^±f±i - n) ^ ^ ^ 

which again implies /3 = —2 — 1/z/ as in Eq. (3.26) for large(r) times, when the quantum correction 
starts to become important (since the left hand side of Einstein's equation always goes like l/t^, 
no matter what the value for a is, at least for k=0). 
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4 Covariant d'Alembertian on Tensor Functions 



Next we will examine the full effective field equations (as opposed to just their trace part) of 
Eq. (2.13) with A = 0, 

R^,v -Ig^^uR = SirG {l + A{a)) T^, (4.1) 
Here the d'Alembertian operator 

□ = g^^-V^V, (4.2) 

acts on a second rank tensor, 

and would thus seem to require the calculation of 1920 terms, of which fortunately many vanish by 
symmetry. Next assume that T^^ has the perfect fluid form, for which one obtains 

(□T^.),, = 6 w) + p{t)\ - 3 Pit) III - m 

(n^^-U = r'{^-kr')sm.'e{nT^,X^ (4.4) 

with the remaining components equal to zero. Note that a non-vanishing pressure contribution is 
generated in the effective field equations, even if one assumes initially a pressureless fluid, p{t) = 0. 
As before, repeated applications of the d'Alembertian □ to the above expressions leads to rapidly 
escalating complexity (for example, eighteen distinct terms are generated by for each of the above 
contributions), which can only be tamed by introducing some further simplifying assumptions. In 
the following we will therefore assume that T^;^ has the perfect fluid form appropriate for non- 
relativistic matter, with a power law behavior for the density, p{t) = po t^, and p{t) = 0. Thus all 
components of T^i^ vanish in the fluid's rest frame, except the tt one, which is simply p{t). Setting 
k = and R{t) = Rq one then finds 

(□ T^v)u = (6 - _ 3 a /? + /?) Po t^~'' 

{^T^u\r = 2I^t^'^o? p^t^-'' (4.5) 
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which again shows that the tt and rr components get mixed by the action of the □ operator, and 
that a non-vanishing rr component gets generated, even though it was not originahy present. 

Higher powers of the d'Alembertian □ acting on T^i, can then be computed as well, but it 
is easier to introduce the slightly more general auxiliary diagonal tensor Vf^u with components 
Vtt = Pot^, Vrr = Pif^j Vee = Vrr ^nd Vipip = r^ sin^ 9 Vrr, with 7 an arbitrary power. One then 
finds 

(□ V^,)^^ = 2Rl t^'^a^ po tf^-^ + (4 + a (7 - 2) - 7 (7 - 1)) pi f^'^ (4.6) 

as well as 

= r''sm''e{DV^,X^ (4.7) 

and zero for the remaining components. The above expressions can then be used conveniently to 
generate acting on T^^, to any desired power n. But since the problem at each step involves a 
two by two matrix acting on the energy momentum tensor, it would seem rather complicated to 
get a closed form solution for arbitrary n. But a comparison with the left hand (gravitational) side 
of the effective field equation, which always behaves like ~ for A; = 0, shows that in fact a 
solution can only be achieved at order provided the exponent /? satisfies (3 = —2 + 2n, or since 
n = -l/{2u), 

P = -2 - 1/u (4.8) 

as was found previously from the trace equation, Eqs. (2.19) and (3.26). As a result one obtains a 
much simpler set of expressions, which now read 

^ Ga^pot-^ (4.9) 

(tfr^^)^^ ^ 12a2(a- l)(4a + l)pot"^ (4.10) 
(c^T^u)^^ 48a2(a-l)(4a + l)(2a2-3a-3)por^ (4.11) 
(p^Ti^i.)^^ 96a2(a- l)(4a + l)(2a2 -3a-3)(4a2 _9a_ i5)pQt-2 (4.12) 
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(n^T^^^^^ 768a2(a-l)(4a + l)(2a2-3a-3)(4a2_9a-15)(a2-3a-7)/9ot~^ (4.13) 

etc., here for powers = 1 to n = 5 respectively, and with /3 changing with n in accordance with 
Eq. (4.8). For general n one can then write 

(□"T^,)^^ ^ Cttia,u)pot-^ (4.14) 

and similarly for the rr component 

{a^T^.\r ^ Crr(,a,u)Rlt^''pot-^ (4.15) 

But remarkably (see also Eq. (4.4) ) one finds for the two coefficients the simple identity 

Crr{a,iy) = lctt{a,u) (4.16) 

as well as c^g = Crr and c^^ = sin^ 9crr- Then for large times, when the quantum correction 
starts to become important, the tt and rr field equations reduce to 

3a2t-2 = 8TrGcttia,u)pot-^ (4.17) 

and 

- a{3a - 2)Rlf°'-'^ = 8nG Crr{a,i^) Rlt^"" pot-"^ (4.18) 

respectively. But the identity Crr = ^ cu implies, simply from the consistency of the tt and rr 
effective field equations at large times, 

Crrioi, l^) _ I 3q; — 2 



ctt{a,v) ^ 3q; 
whose only possible solution finally gives the second sought-for result, namely 



(4.19) 



a = \ (4.20) 

For the specific value of a = | one can then show that the coefficients cu obey the recursion relation 

(c«)n = -(4n2 - 7n+l)(c«)„_i (4.21) 

with initial condition [ctt)n=i = 3/2. Consequently a closed form expression for ctt and Crr = Ctt/'i 
can be written down, either in terms of the Pochhammer symbol (x)„ = x{x + 1) . . . {x + n — 1) = 
V{x + n)/V{x), or more directly in terms of ratios of Gamma functions as 



ctt{a = \,n = -l/2v) = 3(-l)"+^2-=^+^" ^ ^ 8^ 'V ^^ 22) 
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Still, the above expression does not seem to be particularly illuminating at this point, except for 
providing an explicit proof that the coefficients cu and Crr exist and are finite for specific values of 
n, such as n = —l/2v = —3/2. 

One might worry at this point whether the above solution is consistent with covariant energy 
conservation. With the assumed form for T^j^ it is easy to check that energy conservation yields 
for the t component 

(V^ (□"T^.)), ^ - ((3a + /? + l/v)cu + 2,acrr) pot^+'/"~' = (4.23) 

when evaluated for n = —l/2v, and zero for the remaining three spatial components. But from 
the solution for the matter density p{t) at large times one has /? = — 2 — so the above zero 
condition gives again Crr/c-tt = —(3a — 2)/3a, exactly the same relationship previously implied by 
the consistency of the tt and rr field equations. 

In conclusions the values for a = 1/2 of Eq. (4.20) and (3 = —2 — 1/v oi Eq. (4.8), determined 
from the effective field equations at large times, are found to be consistent with both the field 
equations and covariant energy conservation. More importantly, the above solution is also consistent 
with what was found previously by looking at the trace of the effective field equations, Eq. (2.19), 
which also implied the result /? = — 2 — 1/z/, Eq. (3.26). 

Together these results imply that for sufficiently large times the scale factor R{t) behaves as 

R{t) ~ r ~ (4,24) 

and the density p{t) as 

p{t) ~ ~ r^-V- ^ (i?(t))-2(2+V-) (4.25) 

Thus the density decreases significantly faster in time than the classical value (p(t) ~ i~^), again 
a signature of an accelerating expansion at later times. 

It is amusing to note that the vacuum polarization term we have been discussing so far behaves 
very much like a positive pressure term, as should already have been clear from the fact that 
the covariant d'Alembertian g^'^'V causes, in the RW metric case, a mixing of the tt and rr 
components in the field equations. Furthermore, within the classical FRW model, the value a = 1/2 
corresponds to an equation-of-state parameter a; = 1/3 in Eq. (A.20), with 

2 

a = — (4.26) 

3(1 +co) ^ ' 

where p{t) = lo p{t), and which is therefore characteristic of radiation. Thus one can visualize the 

covariant gravitational vacuum polarization contribution as behaving to some extent like classical 

radiation, here in the form of a dilute gas of virtual gravitons. 
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5 Conclusions 



The main results of this paper are the effective field equations of Eq. (2.13), 

R^^v -\9^^uR + ^g„v = SttG (1 + A{n)) r^, (5.1) 

their trace in (Eq. (2.19)), and the solution for the trace and full equations for the specific case of 
the RW metric and A = outlined in Sections 3 and 4, respectively. 

The combined results for the density p{t) ~ Potl^, namely /3 = — 2 — 1/z/for large times 
(Eqs. (3.26) and (4.8)), and for the scale factor R{t) ~ Rot"', namely a = 1/2 (Eq. (4.20)) again 
for large times, imply that for A = and for sufficiently large times the density falls off as 

Pit) ~ r^-V- ^ (i2(t))-2(2+i/-) (5.2) 

Thus the matter density decreases significantly faster in time than predicted by the classical value 
{p{t) ~ t~'^), a signature of an accelerating expansion at later times. 

Within the Friedmann-Robertson- Walker (FRW) framework the gravitational vacuum polar- 
ization term behaves in many ways like a positive pressure term. The value a = 1/2 corresponds 
to w = 1/3 in Eq. (A.20), 

a = — r (5.3) 

3(1 + u) ^ ^ 

where we have taken the pressure and density to be related by p{t) = uj p{t), and which is there- 
fore characteristic of radiation. One can therefore visualize the gravitational vacuum polarization 
contribution as behaving like ordinary radiation, in the form of a dilute virtual graviton gas. It 
should be emphasized though that the relationship between density p{t) and scale factor R{t) is 
very different from the classical case. 

The results of Section 4 show that the effective Priedmann equations for a universe filled with 
non-relativistic matter (p=0) have the following appearance 



R^t) ' \R{t) J 3 ' 3 

SttG 



1 + it/O'^^ + O (mf^-^)] Pit) + I A (5.4) 



R'^it) \Rit) I Rit) 



c^it/^y-^ + O(it/0'/-^)] pit)+k (5.5) 
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with the running G appropriate for the RW metric appearing exphcitly in the first equation, ^ 



G{t) = G[i + c^(t/^)V- + o((i/^)2/.) 



(5.6) 



and used, in the second equation, the resuh cu = 3c,.r of Eq. (4.16). We have also restored 
the cosmological constant term, with a scaled cosmological constant A ~ l/C^- One can therefore 
sensibly talk about an effective density 

G{t) 



Peffit) 

and an effective pressure 



G 



Pit) = l + c^it/O^'^ + 



Pit) 



Peffit) = l(^-^] Pit) 



c^t/O 



1/u 



+ 



Pit) 



(5.7) 



(5.8) 



with. Peffit)/ Peffit) = ^(G(t) — G) /Git). Strictly speaking, the above results can only be proven if 
one assumes that the pressure's time dependence is given by a power law (as discussed in Section 
4). 

Equivalently, substituting t « Q;i?(i)/i?(t), one can, as an example, re- write the first Friedman 
equation as 



k 



' Rit) \ _ SttG 



+ 



(ill' 



+ 



Pit) + u 



(5.9) 



The effective Friedman equations of Eqs. (5.4) and (5.5) also bear a superficial degree of resem- 
blance to what might be obtained in scalar-tensor theories of gravity [25, 26, 27] (for recent reviews 
and further references see [28, 29]), 

1 



levrG 



fi<P)R - \g>'''d^4>du^ - Vi4>) 



+ S, 



matter 



(5.10) 



^Corrections to the above formulae are expected to be fixed by higher order terms in the renormahzation group 
/3-function. In the vicinity of the fixed point at Gc one writes 



p(G) = poiG - Gc) + /3i (G - Gcf + 



and obtains then by integration 



= C^A |exp{- J 



1/0' (Gc) 



with an exponent v given by /3'(Gc) = (3o = —1/v, Cm a numeric constant and A the ultraviolet cutoff. After 
replacing A — » 1/r and G(A) — > G(r) one finds for the scale dependence of G 



G(r) 



1 + 



{l + l3xGcv)Cli" \^ 



2/u 



Note that fio = < Q, and that for /3i < the second correction term is positive as well. If one restricts oneself 

to the lowest order term, valid in the vicinity of the ultraviolet fixed point, then for a given static source of mass M 
one has for the gravitational potential the additional contribution 5V{r) ~ {2MG/^^)r^ for u = 1/3, as discussed 
in [17]. 
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where the gravitational Lagrangian is some arbitrary function of the scalar curvature [30]. It is 
also well known that often these theories can be re-formulated in terms of ordinary Einstein gravity 
coupled appropriately to a scalar field [31]. In the FRW case one has for the scalar curvature in 
terms of the scale factor 

R = 6(k + B?{t) + R{t)R{t)) /R^{t) (5.11) 
and therefore for /c = and R{t) = Rof^, 

R = M2a^ (5^2) 

The quantum correction in Eq. (5.4) is therefore, at this level, indistinguishable from an inverse 
curvature term of the type (^^ i?^)"^/^'". But the resemblance is seen here merely as an artifact 
due to the particularly simple form of the RW metric, with the coincidence of several curvature 
invariants (see for ex. Eqs. (3.8) and (3.12) ) not expected to be true in general. 

Finally let us note that the effective field equations incorporating a vacuum-polarization-driven 
running of G, Eq. (2.13)), could potentially run into serious difficulties with experimental con- 
straints on the time variability of G. These have recently been summarized in [32, 33, 34, 35, 36], 
where it is argued on the basis of detailed studies of the cosmic background anisotropy that the 
variation of G at the recombination epoch is constrained as \G{z = Zrec) — GoI/Gq < 0.05(2cj). So- 
lar system measurements also severely restrict the time variation of Newton's constant to \G/G\ < 
10-i2yr"^ [32]. It would seem though that these constraints can still be accommodated provided 
the scale ^ entering the effective field equations of Eq. (2.13) is chosen to be sufficiently large, at 
least of the order of ^ > 3H~^, given that in the present model one has \G/G\ ~ j; c^t^^'^~^ /^^^^ 
and therefore |G/G| ~ Sc^t'^/^^ for u = l/3. 
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Appendix 



A Classical Field Equations and Conventions 

This appendix is mostly about notation, but also collects a few simple results used extensively 
in the rest of the paper. We will write the Robertson- Walker (RW) metric as 

ds"^ = -dt^ + R^it) 1^3^+^^ (d^^ + sin^^V^l (A.l) 

and note that with this choice of signature (i.e. a minus sign for the dt^ term), □ is a positive 
operator (on functions of t). Also ^/g = \J — det((7) = sin^i?^(t)/-\/l — hr"^. 
The energy momentum tensor for a perfect fluid is 

^ixv = (pit) + p{t))u^u^ + gf,^p{t) (A.2) 

giving in the fluid's rest frame T^i, = diag{p^pB? / (1 — A;r^), r'^pE?, sin^ 6pR^), with trace 

T/ = 3pit)-pit) (A.3) 

The field equations are then written as 

R^.-l9,..R = SttGT^, (A.4) 
The tt component of the Einstein tensor reads 

3(k + R\t)) /R\t) (A.5) 

while the rr component is 



k + if{t) + 2 R{t) R{t)^ ( A.6) 



1 — kr^ 
and the 9 — 6 component 

-r^ (^k + if{t) + 2R{t)R{t)j (A.7) 

and finally the <p — <p component 

-r^ sin^ e(k + R?{t) + 2 R{t) R{t)) (A.8) 

The scalar curvature is simply 

Q(k + R?{t)+R{t)R{t)) /R^{t) (A.9) 
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Thus the tt component of the Einstein equation becomes 

3(k + if{t))/R^{t) = 8TrGp(t) (A.IO) 

while the rr component reads 

■Y^^(k + ifit) + 2Rit)R{t)) = 87rG^-^j^p{t)R'^{t) (A.ll) 
The trace equation is 

6(^k + R^{t) + R{t)R{tfj /R^{t) = SttG {p{t) - 3p{t)) (A.12) 

Covariant conservation of the energy momentum tensor, T^u = imphes a definite relationship 
between R{t), p{t) and p{t), which reads 

m + 3 (p(t) + p{t)) {m/R{t)) = (A.13) 

(and which the tensor of Eq. (A. 2) in its most general form does not satisfy). 

Next consider the case A; = (spatially flat) and p = (non-relativistic matter) . If R{t) = Rq t° 
and p{t) = pat^, then the tt field equation 

^ = MGpot^ (A. 14) 

implies (3 = —2 and = SirGpo/S, while the rr field equation 

- a (3a - 2) i?^ t^"-^ ^ q 

implies a = 2/3. Also both of these together imply 

p(t) ~ ^ (t2/3)-3 _ i/R(tf (A.16) 

The trace equation now reads 

^^^^^ = 8.Gpot^ (A.17) 

and implies again f3 = —2 and 6a (2a — 1) = SvrGpo- The latter combined with the tt equation 
gives 3a^ = 6a (2a — 1), or again a = 2/3. Finally covariant energy conservation implies 

{3a + P)pot^ = (A. 18) 

or 3a + /3 = 0, which does not add to what already comes out of the tt and rr (or, equivalently, 
tt and trace) equations, but is consistent with it. In conclusion the tt and rr (or tt and trace) 
equations are sufficient to determine both a and /3 . 
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The case of non-vanishing pressure can be dealt with in the same way. In most instances one is 
interested in a fairly simple equation of state p{t) = ujp{t), with lo a constant. For non-relativistic 
matter a; = 0, for radiation lo = 1/3, while the cosmological term can be modeled by w = —1. The 
consistency of the tt and rr equations now requires 

a {3a - 2) 



3a2 

which gives 



-LO (A.19) 



a = — (A.20) 

3(1 + Lo) ^ ^ 

for — l<c<;<l/3. Furthermore from the covariant energy conservation law one has 

3{1 + Lo)a + P = (A.21) 
which implies /? = —2 again. Therefore 

R{t) ~ t2/3(i+^) ^(i) ^ (A.22) 

These results are well known and have been collected here for convenient reference. 



B Scale Transformations and Gravitational Functional Integral 

Consider the (Euclidean) Einstein-Hilbert action with a cosmological term 

Ie = ■^oA'^y" dx^ - kqA^ j dxy/gR (B.l) 

Here Aq is the bare cosmological constant, kq = I/IGttGq with Go the bare Newton's constant. 
Also, and in this section only, we follow customary notation used in cutoff field theories and denote 
by A an ultraviolet cutoff, not to be confused with the scaled cosmological constant. The natural 
expectation is for the bare microscopic, dimensionless couplings to have magnitudes of order one 
in units of the cutoff, Aq ~ kq ~ 0(1). Next one can rescale the metric 

g',. = Vx'og.u g'^"" = -^^g^" (B.2) 

to obtain 

Je = A^y dx^ - J dx^R' (B.3) 
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Next consider the fact that the (Euchdean) Feynman path integral 

A2 



^ = J dfJ^ig] exp ^- J dx y/g (^Xo A'^ 



R 



includes a functional integration over all metrics, with functional measure [37, 38] 

dn[g] = 



l[{detG)2l[dg^4x) 



n {g{-)r-'''''^'"' n %^(^ 

X IJ,>V 



D=4 



(B.4) 



n n dgfiuix) 

(B.5) 



with the super-metric over metric deformations given by 

GrM->«/3(5(a;)) = 1 (g{x)f^U-{x)g''>'ix)+g>^^{x)g''''ix) + Xg''''ix)g''fi, 



(B.6) 



For our purposes it will be sufficient to note that under a rescaling of the metric the functional 
measure only picks up an irrelevant multiplicative constant. Such a constant automatically drops 
out when computing averages. Equivalently one can view a rescaling of the metric as simply a 
redefinition of the ultraviolet cutofi" A, A ^ -^y^^- a consequence, the non-trivial part of the 
functional integral over metrics only depends on Aq and kq through the dimensionless combination 
i^q/\/Xq = l/(167rGo\/Ao). The existence of an ultraviolet fixed point is then entirely controlled by 
this dimensionless parameter only, both on the lattice [8, 17] and in the continuum [11, 15]. It is 
the only relevant (as opposed to marginal or irrelevant, in statistical mechanics parlance) scaling 
variable in the pure gravity case, in the sense of having only one positive (growing) eigenvalue of 
the linearized renormalization group transformation in the vicinity of the fixed point. 

The parameter Aq controls the overall scale of the problem (the volume of space-time) , while the 
kq term provides the necessary derivative or coupling term. Since the total volume of space-time 
can hardly be considered a physical observable, quantum averages arc computed by dividing out 
by the total space-time volume. For example, for the quantum expectation value of the Ricci scalar 
one writes 

< J dx y/gixj R{x) > ^g^^ 



< J dx ^g{x) > 

Without any loss of generality one can therefore fix the overall scale in terms of the ultraviolet 
cutoff, and set the bare cosmological constant Aq equal to one in units of the ultraviolet cutoff. 

The addition of matter field prompts one to do some further rescalings. Thus for a scalar field 
with action 

Is = k j dx^ [g''''d^<t>d,(t> + m^</)' + R<t>''} (B.8) 

^"Thoso considerations arc not dissimilar from the case of a sclf-intcracting scalar field where one might want to 
introduce three couplings for the kinetic term, the mass term and the quartic coupling term, respectively. A simple 
rescaling of the field would then reveal that only two coupling ratios are in fact relevant. 
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and functional measure 



(B.9) 



the metric rescaling is to be followed by a field rescaling 



1/4 



(B.IO) 



with the only surviving change being a rescaling of the bare mass mo vn^j \q . The scalar func- 
tional measure acquires an irrelevant multiplicative factor which does not affect quantum averages. 
The bare mass rescaling is of course ineffectual if the fields arc masslcss to begin with. 

The same set of considerations apply as well to the Euclidean lattice [39, 40] regularized version 
of Eq. (B.l), which now reads [41, 42] 



where, as is customary, the lattice ultraviolet cutoff is set equal to one (i.e. all lengths and masses 
are measured in units of the cutoff). It is known that convergence of the Euclidean lattice functional 
integral requires a positive bare cosmological constant Aq > [41, 42, 43]. 

The coupling A should really not be allowed to "run", as the overall space-time volume is 
intended to be fixed, not to be itself rescaled under a renormalization group transformation. Indeed, 
in the spirit of Wilson [4], a renormalization group transformation allows a description of the 
original physical system in terms of a new course grained Hamiltonian, whose new operators are 
interpreted as describing averages of the original system on the finest scale, but within the same 
physical volume. This new effective Hamiltonian is still supposed to describe the original physical 
system, but does so more economically in terms of a reduced set of degrees of freedom. 

The pure gravity theory depends only on one coupling (the dimensionless G), and only that 
coupling is allowed to run. This is also, to some extent, implicit in the correct definition of 
gravitational averages, for example in Eq. (B.7). Physical observable averages such as the one 
in Eq. (B.7) in general have some rather non-trivial dependence on the bare coupling Go, more so 
in the presence of an ultraviolet fixed point. Renormalization in the vicinity of the ultraviolet fixed 
point invariably leads to the introduction of a new dynamically generated, non-perturbative scale 



II = XoY.yh{l'')-2KoY.6h{l^)Ah{l^) 



(B.ll) 



h h 



and 




(B.12) 



for G > Gc 




(B.13) 
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with an exponent related to the derivative of the beta function at the fixed point 

P'iGc) = -Xjv . (B.14) 

The overah size of this new scale ^ is controlled by the distance from the fixed point G — Gc, which 
can be made arbitrarily small (in the Regge lattice theory one finds for the critical coupling, in 
units of the ultraviolet cutoff, Gc ~ 0.626, and for the exponent v k, 0.33) 
Thus a result such as 



< j dx ^/gjx) R{x) > 2/^ ^ ^7I^ a2-7 1 m i 

<!dx^)> ~ ^ ^''-'''^ ~ ^ > ^^-^'^ 

referring here to an average curvature on the largest observable scales (with v and 7 some positive 
exponents) does not presumably allow one to state whether the average curvature is large or small 
at large distances (that would clearly depend on the choice of G — Gc and the cutoff A). But 
it does establish a definite relationship between the fundamental scale ^ in Eq. (B.13) and say 
the scale of the curvature at the largest scales, Eq. (B.15), as well as with any other observable 
involving G — Gc or ^. It is the latter curvature that most likely should be identified with a 
physical, astrophysically measurable, macroscopic cosmological constant (and not in any way with 
Ao). While it is natural to assume for the curvature measured on the largest distance scales (for 
example via the parallel transport of vectors along very large loops) that 1/^^, and therefore 
7 = 2, it has proven difficult so far to establish such a result in the lattice theory, due to the great 
technical difficulties involved in measuring small invariant correlations at large geodesic distances 
[44]. 



C Effective Action Variation 



In this section we will consider the effective gravitational action of Eq. (2.10), 

I = J dx^VR (1 - A{D)) Vr (C.l) 

and compute its variation. One needs the following elementary variations 

S^-VR-il - A{n)) . /r + ^ . 5Vr ■ (1 - A{d)) ■ Vr 
+ ^■VR-d{l-A{n))-VR + ^.VR.{l-A{n)).sVR (C.2) 

^^Pursuing the analogy with Quantum Chromodynamics, wc note that there the non-perturbative gluon condensate 
depends in a nontrivial way on the corresponding confinement scale parameter, as < F^i, ■ F'^" >« (250MeV')* ~ 
with ^Q^o ~ Ajgg-. 
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Using the identity 

SV9 = -hV9g^.u6g'''' (C.3) 

as well as Va g^u = one then has 

- n^VRA{^)^{S^)^/R + ^Vr{1 - A{a))S^/R (C.4) 
Next use is made of the definition of the Ricci scalar, 

SR = g'^'dR^, + R^.dg'^'' (C.5) 
For the variation of the affine connection one has 

ST%, = [V^%, + V,5g^^ - V^%.] (C.6) 

or, equivalently, 

<^rV = -5 [V^C^.A'^^"^) + ^Ag^^^Sg'^^) - Vpig^^g^xg^'^Sg'^'')] (C.7) 
and therefore for the variation of the Ricci tensor 

SR^,. = Va (<5rv) - v^((5r"„j (C.8) 

from which it follows that 

g^^^ SR^, = V^V. ( - 6g^^ + g'^'^ g^p Sg^") = g^p □ Sg^f" - V(^V,) Sg^^ (C.9) 

which is one of the required variations in Eq. (C.4). The second term on the right hand side of the 
last equation is a total derivative the ordinary Einstein case, but it needs to be kept here. Note 
also that in general □ 7^ V^D, and that O g^^,^ = but f^Sg^iy / 0. For the variation of the 
covariant d'Alembertian 

Su = Sg'^'^V^V, - g^^STl^Va (C.IO) 
one needs the variation of F^^ given by Eq. (C.6), which then gives 

6D = Sg^'' V. + Sg^'' V. - ^ g''" g^p Sg'^" V. (C.ll) 

Here (or at the end) one also needs to properly symmetrize the result for the variation of □, 

n 

5{U'^) (C-12) 



fe=i 
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Next several integrations by parts, involving both the operator (with integer n) and well as 
the operator g^u^ — V(^V;y), have to be performed in order to isolate the 6g^'^ term. This follows 
from / -y^g W/^V^ = J yRj [1/ d^j ^JgV^^ = which allows us to repeatedly integrate by parts 
and move some covariant derivatives around. In general one has to be careful about the ordering 
of covariant derivatives, whose commutator is in general non-zero in accord with the Ricci identity 

* j 
with the a index in T in the i-th position in the first term, and in the j-th position in the second 
term. The term involving the variation of the covariant d'Alembertian □ then gives 

-n (v,V.x/i?) [^^) -n {v,VR) {v^^^)+\ng,. (v«v^) g'^^ (^P^^ 

(C.14) ■ 

which again needs to be symmetrized with respect to \/R ^ \/R, in the way described above. 
After adding the remaining terms, the effective field equations become 

- hg^^uR) (l - -^A{u)^^R)J - (g^,u - V(^V,)) 1^-1= A{u)^ 



n 



+ \ng^, (VaVR) ff"^ (v^^/r) = SttGT^, (C.15) 

where again the last three terms need to be properly symmetrized in \fR ^ ^/R, as described 
above. 

Taking the covariant divergence of the l.h.s gives zero for some of the terms, while the remaining 
terms give 

1 



which has to vanish due to the invariance of the original non-local action. (Again the last term 
needs to be symmetrized in D"""^ \fR <->■ s/R). 

The above derivation can be slightly generalized to an action of the form 

I = I dx^R'--{l - A{n))R- (C.17) 
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with a a parameter between zero and one (with the previous case corresponding to a = 1/2). Then 
for the field equations one obtains an expression of the type 

R^v - ^f^uR + \g^uRR-''A{u)R'^ - R^, ((1 - a) 72"" ^(n) i?" + ai?"-^ A(n)i?i-") 
- (g^^uO - V(^V,)) ((i-a)i?-"A(n)i?« + aR^-'A{a)R^-^) 
-n (V^V.i?") (^^R'-^ - n (V^i?") (v,^R^--^ 

+ \ng^, iV^R'^) g-"" (^Vp^i?^"") = SttGT^, (C.18) 

(where again the last term needs to be symmetrized) which shows that the choice of either a = 1 
or Q = is a bit problematic. 

One final question remains, namely what is the relationship between the above effective field 
equations, Eq. (C.15) or Eq. (C.18), and the clearly more economical field equations proposed in 
Eq. (2.13). Obviously the equations obtained above from the variational principle are much more 
complicated. They contain a number of non-trivial terms, some of which are reminiscent of the 
1 + A{B) term, and others with a completely different structure (such as the g^i^ □ — V(juVi^) 
term). It is of course possible that when restricted to specific metrics, such as the RW one, the 
two sets of equations will ultimately give similar results, but in general this remains a largely open 
question. One possibility is that both sets of effective field equations describe the same running 
of the gravitational coupling, up to curvature squared (higher derivative) terms, which become 
irrelevant at very large distances. 
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